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The direct correlation function and the (static) structure factor for a seven-dimensional hard- 
sphere fluid are considered. Analytical results for these quantities are derived within the Percus- 
Yevick theory 



Interest in studying hard-core fluids in different dimen- 
sionalities is tied to the notion that, since these systems 
share some characteristic features such as the existence 
of a freezing transition, by considering a higher dimen- 
sion [in which the mathematics may sometimes become 
simpler (see for instance Ref. [ll)] one may gain insight 
into common phenomenology which is either untractable 
or rather difhcult in two or three dimensions. Therefore 
it is not surprising that many studies have been devoted 
to hard-core fluids in dimensions higher than three, in- 
cluding fairly recent ones, mostly involving virial coeffi- 
cients and/or equations of state (for a non exhaustive but 
hopefully representative list see Refs.Siiili). On the 
other hand, the importance of the Percus- Yevick (PY) 
theory^ for the structural and thermodynamic properties 
of liquids rests on the fact that it is exactly solvable in 
the case of hard-core systems in odd dimensions i^i^ In 
fact, the direct correlation function c(r) and the static 
structure factor S{q) have been derived analytically in 
one, three, and five dimensions (see Ref. d and references 
therein). The major aim of this Note is to provide the 
explicit expressions for the functions c(r) and S{q) of the 
7D hard-sphere fluid in the PY theory. This work ex- 
tends and complements our previous paper— in which we 
concentrated on the thermodynamic properties and the 
virial coefficients of the same system, also within the PY 
approximation, and on the comparison of some of these 
results with our own simulation data. It must be em- 
phasized that the highest dimensionality for which com- 
pletely analytical results may be derived from the PY 
theory is precisely d — 7. Beyond that dimensionality, 
numerical work is required at one stage or another. 

Following Ref. d, the direct correlation function of the 
7D hard-sphere fluid has, in the PY theory, the polyno- 
mial form 



c{r) = (co + cir + c^^r^ + c^r^ + c^r'^) 8(1 - r), (1) 



where the diameter of a sphere is cr = 1 and Q{x) is Heav- 
iside's step function. The coefficients Ci {i = 0, 1, 3, 5, 7) 
depend on the packing fraction r] = (7r^/840)yO (where p 
is the number density) . After some algebra, one finds the 



following structure: 



^ " j=0 



(2) 



Here, Q^^^ is the physical solution of a quartic equation;^ 
which reads 



)(0) 



336077(1 - 77) 



- , (3) 



where 



Ti = ^(3 + 3677 + 10772) + is. 



(4) 



T2 = 2Ti - is - (3 - 11477 - 2227^2 - 1077^) (5) 



5 = ^1/^ + (3 + 367; +10772)^^-1/3, (6) 
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Ai + ^jAf- (3 + 3677 + IO772)" 



(7) 



Ai ^ 459 + 4374O77 - 464477^ + 1652477^ - 35794877^ 
+16092077^ + 23300r;^ (8) 

The coefficients tij appearing in Eq. ^ are listed in Table 
m It can be checked that c{l^) = —g{l^), as expected 
from the continuity of the indirect correlation function 
7(r) = g{r) — c{r) — 1. As an illustration, the top panel 
of Fig. [1] shows c(r) for 77 = 0.01, 0.02, . . . , 0.07. 

Now we turn to the structure factor S{q). According 
to the Ornstein-Zernike equation, S{q) is related to the 
Fourier transform 'c{q) of c(r) by 



S{q) 



1 



1 - pc{q) 
In 7D, the Fourier transform is^ 



/•oo 

c(g) = (2^)^/2g-5/2 / drr"^J^/^{qr)c{r), 
Jo 



(9) 



(10) 
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TABLE I: Expressions for the coefficients t. 



2 toj 

-(1 + 6j;)^ 

1 2688»)(2 + 5)))(1 + 6r;) 

2 - 1053696?)^ (9 + 50?j + 60r;^) 

3 10115481601)^ (3 + 4r))(2 + 5r)) 

4 -1416167424007)'' (3 + 4i))^ 

j *l3 

^1,(4 + 3r))^ 

1 -840t7(4 + 3r)) (1 + 33r; + 15))^) 

2 6720)) (12 + 1212?) + 18783r)^ + 17060r)^ + 3750?;'') 

3 -45158400))^ (3 + 36?) + 10?)^) (l + 33j) + 15j)^) 

4 6322176000?)^ (3 + 36?) + 10?)^)'' 

i t_32 

-f|(20 - 13r)) (12 - 38?) + 33r;'') 

1 -420 (60 - 2759?) + 9163?)^ - 10975?)^ + 3825?; ■*) 

2 3360?) (22041 - 308370?) + 812905?)^ - 882160?)^ + 273950?)*) 

3 -22579200?)^ (3147 - 19924?) + 38685?)^ - 36360?)^ + 9650?)*) 

4 6322176000?)^ (3144 - 125767) + 17541?)^ - 13360?;^ + 2850?)") 

i tsj 

fi (1320 - 5107?) + 6282?)^ - 2446?)^) 

1 252 (330 - 15179?) + 51509?)^ - 59291?)^ + 21945?)*) 

2 -14112?) (17325 - 242482?) + 661205?)^ - 674260?)^ + 226550?)*) 

3 94832640?)^ (2475 - 15664?) + 31953?)^ - 27500?)^ + 8050?)*) 

4 -26553139200?)^ (2475 - 9888?) + 14702?)^ - 10032?)^ + 2400?)*) 

j *7j 

- (20592 - 79789?) + 97872?)^ - 384307)-') 

1 -24 (2574 - 118404?) + 402761?)^ - 4613417)^ + 172695?)*) 

2 47040?) (3861 - 54054?) + 147942?)^ - 149642?)^ + 51060?)*) 

3 -45158400?)^ (3861 - 24453?) + 50151?)^ - 42632?)^ + 12730?)*) 

4 6322176000?)^ (7722 - 30888?) + 46269?)^ - 31056?)^ + 7610?)*) 



TABLE ff: Expressions for tfie functions Ci(q). 



2 

q (q' - 15) 

1 g* - 249^ + 48 

3 - 48<)* + 576g^ - 1152 

5 q^ - SOq" + 2400ij* - 288009^ + 57600 

7 f/'" - 120()** + 6720()'' - 201600()* + 2419200g^ - 4838400 



where J5/2{z) = ^2/7rz~^/^ [(z^ — 3) sin z — 3z cos z] is 
a Bessel function of the first kind. Inserting Eq. ^ into 
Eq. (jlO[) and performing the integrals one gets 

c{q) = co^Aolg) + cit/jiiq) + csipsiq) + C5'05(g) + criJTiq), 

(11) 

where 

V'7(g) = 1677 V^^+'^ [C^iq) COSq-S,iq)smq-C,{0)]. 

(12) 

The functions Ci{q) are hsted in Tabic HH while the func- 
tions Si{q) are 

5,(g) = -fc.(q) + 3<z2+\ (13) 
aq 

Insertion of Eq. pT|) into Eq. gives the structure fac- 
tor. This function is plotted in the bottom panel of Fig.[T] 



In view of the rather satisfactory agreement that has 

been observed between the analytical PY results for the 

compressibility factor of the 7D hard-sphere fluid and 

the simulation data,^'^ one wonders whether a similar 
0| • , • , • , . I . I 
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FIG. 1: Plot of the direct correlation function c(r) (top panel) 
and of the structure factor S{q) (bottom panel) of the 7D fluid 
of hard hyperspheres in the PY theory. The packing fractions 
corresponding to the curves are, from top to bottom in the 
left end, 77 = 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, and 0.07. 



agreement will hold for the structural properties. That 
this is indeed the case has been recently confirmed by 
Bishop and his coworkers;^ who have found a very good 
agreement between their simulation results and the ones 
computed using our analytical formulae. 
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